Jamming in Systems Composed of Frictionless Ellipse-Shaped Particles 
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We study the structural and mechanical properties of jammed ellipse packings, and find that the nature of the 
jamming transition in these systems is fundamentally different from that for spherical particles. Ellipse packings 
are generically hypostatic with more degrees of freedom than constraints. The spectra of low energy excitations 
possess two gaps and three distinct branches over a range of aspect ratios. In the zero compression limit, the 
energy of the modes in the lowest branch increases quartically with deformation amplitude, and the density of 
states possesses a 5-function at zero frequency. We identify scaling relations that collapse the low-frequency 
part of the spectra for different aspect ratios. Finally, we find that the degree of hypostaticity is determined by 
the number of quartic modes of the packing. 

PACS numbers: 61.43.-j, 81.05.Kf, 63.50.Lm, 83.80.Fg 
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A decade ago, Liu and Nagel proposed that jamming transi- 
tions in glassy, granular, and other nonequilibrium amorphous 
systems can be described by the jamming phase diagram iH, 
and that a 'fixed point' (Point J) in the jamming phase diagram 
controls slow dynamics in these systems even far from Point 
J iQl- For model disordered systems composed of friction- 
less, spherical grains interacting via purely repulsive, short- 
range potentials, the packings at Point J are isostatic ^ 0], 
where the number of degrees of freedom exactly matches the 
number of constraints. It has been shown that isostatic pack- 
ings of spherical grains have interesting mechanical proper- 
ties; for example, they possess an abundance of spatially ex- 
tended "floppy modes" of excitation jsl S 01 and non-elastic 
response to applied deformations |@1 . 

However, there have been relatively few theoretical stud- 
ies of jamming in systems with aspherical particles, despite 
the fact that most physical particulate media have grains with 
anisotropic shapes. The introduction of aspherical particles in 
equilibrium systems gives rise to completely new phases of 
matter and critical behavior as evidenced in the field of liquid 
crystals. What is the nature of jamming transitions in nonequi- 
librium systems when the grains are aspherical? For example, 
is there still a special point J in the jamming phase diagram, 
where packings are isostatic, that controls slow dynamics? In 
this letter, we begin to address these important open questions 
by investigating the structural and mechanical properties of 
static packings of frictionless anisotropic particles. 

For a static packing of grains in d spatial dimensions, 
with d — \ rotational and d translational degrees of free- 
dom per grain, total force and torque balance on each grain 
can be satisfied only if the total number of contacts satis- 
fies Nc > Ni = N {2d — I) in the large system limit. Isostatic 
packings satisfy A^^ = Nj, while hypostatic packings possess 
Nc < Nj. In contrast to spherical particle packings, static pack- 
ings of ellipsoidal particles, studied previously by Donev, et 
al. 10], are generically hypostatic and possess higher densi- 
ties without translational and orientational order Experiments 
on packings of ellipsoidal M&M candies have also verified 



these results lfiol] . These previous studies raise several fun- 
damental questions about static packings of ellipsoidal par- 
ticles; for example, why are they hypostatic and what is the 
nature of their low-energy excitations? In this letter, we study 
the low-energy, vibrational excitations of 2D static packings 
of ellipses using a numerical packing-generation algorithm in 
which soft ellipses interact via purely repulsive potentials at 
zero temperature. Our analysis demonstrates the existence of 
two gaps in the vibrational spectrum over a range of aspect 
ratios. The energy of the modes below the first gap increases 
quartically with deformation amplitude along the soft direc- 
tions, and the number of these quartic modes determines the 
degree of hypostaticity of the packings. 

Compression packing-generation protocol We generated 
an extensive set of static packings of ellipses over a range of 
system sizes from = 120 to 1920, in which particles are 
'just touching', using a numerical packing-generation proto- 
col similar to that employed to create static packings of spher- 
ical particles 11211 . We will refer to this protocol as the 
'compression method'. In this method, soft, purely repulsive 
ellipses are first randomly deposited in a square cell with peri- 
odic boundary conditions at a low packing fraction ((|) = 0.5). 
The configurations are successively compressed in small steps 
(A(|) = 10^^) and then relaxed using conjugate energy min- 
imization after each step. Near the jamming packing frac- 
tion (|)y, where the particles are just touching, the configura- 
tions are expanded or compressed by decreasing amounts un- 
til the system has vanishingly small total potential energy per 
particle Vtoi <V < 2ytoi. Vtoi = 10^'^ for most simulations, 
V — L;>; ^(''i;) Summed over all ellipse pairs. 
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and rij is the separation between the centers of mass of ellipses 
/ and j. The separations and orientations m, and uj of the long 
axes of ellipses ; and j determine the Perram and Wertheim 
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FIG. 1: (Color online) Normal mode frequencies CO; from the dy- 
namical matrix (Eq. |3j vs. index sorted by increasing frequency 
for A' = 120 ellipse packings at nine aspect ratios, a = 1.02 (black), 
1.04 (red), 1.06 (green), 1.08 (blue), 1.1 (yellow), 1.2 (brown), 1.4 
(gray), 1.6 (violet), and 1.80 (turquoise). The sorted frequency spec- 
trum possesses three distinct branches (numbered 1, 2, and 3). In the 
inset, we show the scaled frequency Q),/(0* vs. i/i*, which collapses 
the low frequency part of the spectra at 1 (vertical dashed line). 
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where Oq, p, and % depend on A, and the major (minor) axis. 



fl/ (bi), of the ;th ellipse llTT . To determine a,y for each pair, 
minimization of the parameter < A- < 1 must be performed. 
We simulate bidisperse mixtures to inhibit translational and 
orientational order: one-third of the particles are large with 
the major axis 1.4 times that of the small particles ||9[]. Us- 
ing this procedure, we generated an ensemble of at least 100 
ellipse packings, each characterized by the jamming packing 
fraction over a range of aspect ratios from a = 1 to 2. 
We calculated t he g lobal bond-orientational and nematic or- 
der parameters [ISfl, and found no significant ordering with 
order parameters ~ 1 / V^V for all a. 

Vibrational Spectra To determine mechanical properties 
of ellipse packings, we calculate the dynamical matrix 



M,„„ = 



(3) 



where = {x,n,ym,cimQm}, x,,, and y,„ are center of mass 
coordinates of ellipse m, 9m is the angle between u,„ and 
the X-axis, and m,n = 1,- -,A? d. WhenEq.E] is evalu- 
ated for an ellipse packing and diagonalized using periodic 
boundary conditions, in principle one obtains {2d — 1 )A^' — d 
nontrivial vibrational eigenmodes, where N' = N — N,- and 
Nr is the number of 'rattlers' with fewer than d + 1 con- 
tacts. If we assume that all ellipses have the same mass. 
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FIG. 2: (Color online) Change in potential energy AV,- vs. displace- 
ment S along ei for N = 120 and a = 1.5 (/* = 22). AV,- (solid blue) 
for ; = 115 > i* is quadratic in 6. In contrast, for i < i*, AV/ oc for 
5 < 8c, but oc 6"* for 5 > 5c. For / = 24, we show that Sc decreases 
from approximately lO^-' to 10^"* as V,oi varies from 10^'^ (green) 
to 10^'* (purple). The dashed black (red) line has slope two (four). 
Inset: The translational contribution 7] to the sum of the squares of 
the amplitudes of each eigenvector e/ of the dynamical matrix for as- 
pect ratio a = 1.01 (black), 1.20 (red), 1.50 (green), and 2.00 (blue). 

the square roots of the eigenvalues of the dynamical ma- 
trix give the normal mode frequencies co, indexed by ;. We 
denote the normalized eigenvector corresponding to co, by 



ei = Wxi ''^e; • • • '4 ™ '4 " '4 " } with the con- 
straint that e? = 1 . Below, we show the relative contributions 
of translational and orientational components of the eigenvec- 
tors, for example, the translational contribution from mode / 



is T: 
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(e;.)2}and7?, = l-7;-. 



Over a range of aspect ratios, the spectrum CO,, sorted in 
order of increasing frequency, possesses three distinct regimes 
(c/Fig. [T]): (1) modes with indexes i — 2< i*{a) below the 
low-frequency gap, (2) modes with index /* < i — 2 < it = 
{d — l)N, where for a < ttr, there is a second gap at index it, 
and (3) modes with index / — 2 > (We are explicitly not 
including the two trivial modes corresponding to translational 
invariance.) In the inset to Fig. [T] we show that we are able 
to choose aspect ratio dependent scaling factors co* and ;* that 
collapse the low-frequency part of the spectra including the 
first gap. We find that CO* is roughly linear with a — 1, while 
i* possesses two different scaling regimes for a — 1 ^ 1 and 
a > 1 . As demonstrated in the inset to Fig.|2] modes in regions 
(1) and (2) are mainly rotational, whereas those in region (3) 
correspond to mainly translational motion. 

We find that our ellipse packings at finite, but small over- 
compression possess {2d — l)N — d nonzero, positive eigen- 
values 120(1. To rationalize this result with the fact that hard 
ellipse packings are generically hypostatic |@], we perturbed 
our ellipse packings along each of the eigendirections of the 
dynamical matrix over a range of overcompression. If ^0 char- 
acterizes the centers of mass and orientations of the original 
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FIG. 3: Density of vibrational modes O(£0) for ellipse packings with 
A' = 120 for a = 1.0 (black), 1.1 (red), 1.2 (green), 1.4 (blue), 1.7 
(violet), and 2.0 (orange). The low CO peak corresponds to rotational 
modes in frequency regime 2 discussed in the text. The inset magni- 
fies the low CO region for ot near 1: a = 1.001 (thin solid), 1.005 (thin 
dashed), 1.01 (thin dotted), 1.02 (thick sohd), 1.04 (thick dashed), 
and 1.08 (filled circles). The low-frequency peak is sharp and Z?(co) 
possesses a gap in CO for a < a, ~ 1.2. However, the peak broadens 
and connects to Z)(co) at large CO without going to zero for a> a,. 

static ellipse packing, the perturbed configuration obtained af- 
ter a shift by 5 along eigenmode / and relaxation to the nearest 
local energy minimum is ^, = + 5e;. In Fig. |2] we plot the 
change in potential energy, AV/ = V{t,i) — V{t,o), of ellipse 
packings with = 120 at a = 1 .5 arising from a perturbation 
along mode ; as a function of amplitude 8 and for two val- 
ues of overcompression Vtoi. As shown in Fig. |2] for modes 
with indexes in regions (2) and (3) of the frequency spectrum 
AV,- oc for all 8 independent of Vtoi. In contrast, there is 
range 8 > 8c over which modes in region (1) display quartic 
dependence on 8, AV; 8'*, whereas AV, °^ ?r for 8 < 8c . Since 
8c ~ V^J^^ for modes in region (1), quartic behavior will per- 
sist over the entire range of 8 in the zero compression limit. 
Thus, 'just-touching' ellipse packings are stabilized by quar- 
tic rather than quadratic terms in the expansion of the total 
potential energy around the reference packing 

The density of vibrational modes, D(co), obtained from the 
spectrum, co,, shown in Fig. [3] exhibits several key differences 
from that for disk packings First, the plateau at small 
CO characteristic of nearly-isostatic disk packings is replaced 
by a narrow peak at low frequencies. This peak is composed 
of modes in region (2) of the spectrum that display collec- 
tive, primarily rotational motions. For a < ttr (with a? w 1.2 
for = 120), this peak is clearly separated by a gap from 
the broad, high-frequency regime. Note that on a logarithmic 
scale one would also see a peak that corresponds to modes in 
region (1), and is separated by a gap from the modes in region 
(2). As discussed before, this peak shifts to lower co, and nar- 
rows to a 8-function at zero frequency in the zero-compression 
limit. Second, for large aspect ratios (a > 1.5 for = 120), 
D(co) has significantly more structure than that for disk pack- 
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FIG. 4: (Color online) Contact number (z) vs. a for compression 
(triangles) and annealing (Aa = 0.005 (open squares) and 0.05 (filled 
circles)) methods for A' = 480. (z) for annealing diverges from that 
for compression for a > 1.1. The filled diamonds represent (z) = 
6 — i*(a)/{3N) for the compression method data. The inset shows 
the system size dependence of A'/r(a) — A'/(A' + 1) for A' = 120 
(circles), 240 (squares), 960 (upward triangles), and 1920 (downward 
triangles) for the compression data. The red (black) dashed line has 
slope 0.5(1.0). The solid line is an interpolation between power-laws 
of 0.5 at small a— 1 and I.O at large a— 1, which was used to fit (z) 
in the main plot (solid blue line). 

ings for co > 10. These new features will be investigated in 
more detail in future work 12 lil . We have verified that the vi- 
brational properties described here also hold for packings in 
which ellipses interact via the Gay-Berne potential ll22ll with 
a simpler form for the overlap function than in Eq.|2l 

Hypostaticity Our analysis of the structural and mechan- 
ical properties of ellipse packings has led to several novel 
observations: (1) The quartic modes that exist for hard el- 
lipse packings develop a quadratic contribution at finite over- 
compression, and therefore compression stabilizes ellipse 
packings!^; (2) The contribution of these quartic modes to 
the density of vibrational modes decreases as the aspect ra- 
tio increases; and (3) Ellipse packings are generically hypo- 
static, as shown in Fig.|4] with smaller average contact number 
(z) than that predicted from isostatic counting arguments, i.e. 
{z)iso = 6 for ellipses with a > 1. Hypostatic packings have 
fewer contacts than required to satisfy all force and torque bal- 
ance conditions, and thus in some directions of configuration 
space these packings are only quartically, not quadratically 
stable. We find that quartic modes represent collective, pri- 
marily rotational motions of the ellipses, which do not lead to 
cage breaking and particle rearrangements. Thus, we expect 
that if the isostatic counting argument is reformulated so that 
the quartic modes are not constrained, (z) will correspond to 
the minimum number of contacts necessary to constrain the 
quadratic modes. Thus, ellipse packings are isostatic with re- 
spect to only the quadratic modes. 

When each internal degree of freedom in an ellipse packing 
is stabilized, the isostatic conjecture gives A^(z) /2 — 3N— 1. If 
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it is not necessary to constrain the quartic modes, this equation 
can be rewritten as 



N{z 



= 3N-i*{a), 



(4) 



where A^4 = /* (a) — 1 is the number of quartic modes in region 
1 of the frequency spectrum. Thus, by measuring the number 
of quartic modes, we are able to determine (z) (a) as shown 
in Fig.m In the inset we show that A^/;*(a) -A^/(A^+ 1) has 
two power-law regimes: scaling as \/a— 1 (a — 1) for small 
(large) a — 1. In these limits, we obtain 



(5) 



where n = 0.5 (1) for small (large) a — 1 and Ao.s and Ai are 
positive constants. In Fig. IH we use a function that interpo- 
lates between these power-laws and allows us to fit (z) for the 
compression method over the entire range of a. These argu- 
ments imply that (z) = (z)iso in the a ^ °° limit, yet this is 
still an open question. Previous studies have predicted that 
(z) (a) — (z) ( 1 ) scale as Va^~T based on the behavior of the 
pair distribution function g{r) near contact for spherical par- 
ticle packings @]. In contrast, our numerical results demon- 
strate hypostaticity in ellipse packings originates from the ex- 
istence of quartic modes. 

Annealed packings The ellipse packings discussed up to 
this point were generated using the compression method at 
fixed aspect ratio. Using this method, we obtained similar 
(z) (oc) to that obtained previously for hard ellipse packings 
191]. Since ellipse packings are hypostatic, it is in princi- 
ple possible to obtain packings with higher (z) than found in 
Fig. m without increasing the translational or orientational or- 
der To test this, we developed an 'annealing' method, which 
creates static packings by incrementally increasing the aspect 
ratio from a = 1 . We initially create bidisperse disk packings. 
Each disk is then assigned the same aspect ratio 1 + Aa with 
the direction of the long axis chosen randomly. A new ellipse 
packing is formed from this initial state using the compression 
method described above. The ellipses of the new packing are 
elongated again along their defined major axes, and the pro- 
tocol is repeated until a packing with the desired aspect ratio 
is reached. Using this method, ellipse packings can be gener- 
ated that are denser and possess contact numbers much closer 
to Ziso = 6 as shown in Fig.|4] The annealed packings still ex- 
hibit a low-frequency gap, and the location, /*, of this gap can 
be used to predict (z) through Eq. ID The variation of ;* with 
aspect ratio, however, differs from the 'compressed' packings. 
The annealed packings tend to have higher (z), fewer quartic 
modes, and a plateau in (z) at large a. 

In conclusion, we find that aspherical grains qualitatively 
change the nature of jamming at point J. Model systems 
with spherical grains appear to be the exceptional case: they 
are isostatic, and all nontrivial modes of excitation increase 
quadratically with deformation amplitude. In contrast, ellipse 
packings, which are more relevant for real granular media, 
possess quartic modes characterized by collective rotational 



motions. Thus, Landau theories for ellipse packings have a 
4th order term as the lowest nonvanishing contribution to the 
generalized free energy. These results will likely stimulate 
further work on statistical descriptions of packings of aspher- 
ical particles. 
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